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In the decays of beauty to charm hadrons the fun- 
damental force underlying the decay of a 6 quark to a 
c quark is the weak interaction. However, the heavy b 
quark is surrounded by a cloud of light quarks and gluons 
so the strong interaction must be taken into account. In 
the limit of an infinite mass of the b quark, the heavy 
quark decouples from the light degrees of freedom. For 
a finite the decay rates can be computed as a series 
expanded in the small parameter Aqcd l m b, where m& is 
the mass of the b quark and Aqcd is the energy scale of 
the QCD interactions within the hadron. This is known 
as the heavy-quark expansion (HQE) 1]. The applica- 
tion of HQE to the decays of the A° baryon (udb) and 
the beauty mesons (B°,bd; B + ,bu) does not result in an 
identical series. For example, in the (Aqcd /mb) 3 term 
IF-boson exchange contributions are quite different Q, 
leading to a prediction of t(A®)/t(B°) ^ 1. Experimen- 
tal studies of beauty hadron lifetimes therefore help us 
to test the theoretical understanding of the HQE series, 
and consequently the underlying QCD physics. 

Over the past five years theoretical predictions of the 
lifetime ratio r(A°)/r(B°) have not agreed with exper- 
imental values. In 2004 an HQE calculation including 
0(l/m$) effects resulted in r(Ag)/r(S°) = 0.86±0.05 @|. 
This was in good agreement with the 2006 experimental 
world average of 0.804 ± 0.049 0]. In 2006 the CDF 
collaboration reported a measurement @ of the A° life- 
time in the Ag — » J/ipA° channel such that t(A°)/t(B°) 
differed by +2a from the 2006 world average [4J, was sig- 
nificantly higher than the 2004 HQE calculation 0], but 
was compatible with earlier HQE predictions Q . A more 
recent measurement by the D0 collaboration 0] in the 
same channel leads to a value of t(A®)/t(B°) which is 
compatible with both the 2006 world average [4|] and the 
CDF value 0. 

In this paper we present the first measurement of the 
A£ lifetime in a fully hadronic final state. The data sam- 
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pie is produced in pp collisions at y/s = 1.96 TeV at the 
Tevatron and corresponds to an integrated luminosity of 
1.1 fb _1 . We reconstruct A° in the A° — > A+7r~ decay 
channel where A+ subsequently decays as A+ — > pK~ir + . 
Throughout the paper, reference to a specific charge state 
also implies the charge conjugate state. 

The components of the CDF II detector || most rel- 
evant for this analysis are the tracking system and the 
displaced vertex trigger system. The tracking system lies 
within a uniform axial magnetic field of 1.4 T. The inner 
tracking volume is instrumented with either 6 or 7 layers 
of double-sided silicon microstrip detectors up to a radius 
of 28 cm from the beamline 0] . These surround a layer of 
single-sided silicon mounted directly on the beam pipe at 
a radius of 1.5 cm [l(j. This system provides an excellent 
resolution (about 40 /im) on the impact parameter (do), 
which is defined as the distance of closest approach of the 
charged particle to the pp interaction point in the plane 
transverse to the beam direction. The do resolution of 
40 /im includes an approximate 30 /im contribution from 
the uncertainty of the interaction point in the transverse 
plane (added in quadrature). The outer tracking volume 
contains an open-cell drift chamber (COT) up to the ra- 
dius of 137 cm[n|]. 

CDF II employs a three-level trigger system. The ex- 
tremely fast tracker (XFT) [l2| at the first level groups 
COT hits into tracks in the transverse plane. At the 
second level, the silicon vertex trigger (SVT) [l3| adds 
silicon hits to the tracks found by the XFT, improv- 
ing the resolution of the track position and thus allow- 
ing selection based on the transverse displacement from 
the beam line that is measured in real time. The dis- 
placed vertex trigger 14| requires two charged parti- 



cles with momentum transverse to the beam direction 
(pt ) greater than 2 GeV/c, and with impact parame- 
ters in the range 0.12 < \do\ < 1 mm. The intersection 
point of the two particle trajectories must have a trans- 
verse displacement (L xy ) from the interaction point of 
at least 200 p,m. The pair must also have a scalar sum 
Pt(1) + Pt(2) > 5.5 GeV/c. This trigger configuration 
based on a pair of tracks is called the two-track trigger 
(TTT) and is the basis for the collection of many fully 
hadronic bottom and charm decays at CDF. 

We reconstruct a A" candidate via its decay to A+ n~ , 
where the A+ further decays to a pK~ir + final state. All 
four tracks are required to have a sufficient number of 
hits in the tracking detectors for high-quality position 
measurement. Several requirements are imposed to sup- 
press background in the reconstructed sample which are 
optimized using simulated signal and data background 
samples [15(. Each particle must have |g?o| < 1000 fj,m. 
We construct A+ candidates by combining three tracks 
assuming the (pK~ir + ) hypothesis. Thep candidate and 
the ir~ are required to havepx > 2.0 GeV/c. The proton 
Pt must exceed the pt of the ir + from the A+ , which 
has the same charge. This prevents the same pair of 
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tracks being considered both as (p, ir + ) and as (n + ,p). 
The three tracks from the A+ candidate are first con- 
strained to a common vertex in a kinematic fit. Next 
we add a track and construct A° candidates through a 
further kinematic fit, which intersects the fourth track 
with the A+ candidate trajectory. The mass of the A+ 
candidate is constrained to the world average A+ mass 
(2.286 GeV/c 2 ) [lf|. This second kinematic fit allows 
us to calculate ct— L xy cM/pT and its uncertainty, tr c t, 
where c is the speed of light, and t and M are the proper- 
decay time and measured mass of the A° , respectively. 

We apply additional selection requirements in order 
to suppress background. The requirements on ct(A^) > 
250 jitm, its significance ci(A°)/er ct > 10, and |rfo(A°)| < 
80 /zm primarily suppress the background arising from 
random combinations of tracks, many of which originate 
from the primary interaction point (combinatorial back- 
ground). Another important source of background is 
the decay of B mesons with misidentified decay prod- 
ucts. Decays like B° — ► D + ir~ are especially insidi- 
ous since they are abundant (compared to A° — > A+7r~ 
decays) and D + — > K~ii + ii + decays can easily mimic 
the A+ — > K~pir + signature. These backgrounds are 
suppressed by selecting a narrow region of the invari- 
ant mass spectrum of the A+ — ► pK~ir + candidate: 
\m(pK-Tr+) -m(A+) PDG \ < 16 MeV/c 2 . Further D+ 
candidates are removed by a requirement on the ct of A+ 
candidates with respect to the A° vertex, since the A+ 
candidates are usually much shorter lived than the D + 
candidates. We require —70 < ct(A+w.r.t.A°) < 200 /im. 
Lastly, the TTT criteria are confirmed using the recon- 
structed candidate tracks. 

The lifetime of the A° baryon is determined from two 
sequential maximum likelihood fits. The first is a fit to 
the invariant mass of A+7r~ candidates and is used to 
establish the composition of the sample. This gives the 
normalization of each of the fit components for both the 
whole domain of 4.82 < m(Ag) < 7.0 GeV/c 2 , as well 
as the signal region (5.565 < m(Ag) < 5.670 GeV/c 2 ). 
The second fit is an unbinned maximum likelihood fit of 
ctand (T c tin the signal region to extract the A" lifetime 
with the normalizations of each component fixed. 

The invariant mass distribution of A+7r~ candidates is 
shown in Fig. Q] with the fit projection overlaid. Small 
deviations of the model from data below the A° mass do 
not affect the lifetime as they occur outside the signal 
region. The A+7r~ mass distribution is described by sev- 
eral components: the A° — > A+7T - signal, a combinatorial 
background, partially and fully reconstructed B mesons 
that pass the A+7r~ selection criteria, partially recon- 
structed A° decays, and fully reconstructed A° decays 
other than A+7T - (e.g. A° — > A+AT - ). The combinatorial 
background is modeled with an exponentially decreasing 
function of A+7r~ mass. All other components are rep- 
resented in the fit by fixed shapes derived from Monte 
Carlo (MC) simulations [l7| whose relative contributions 



are constrained using data when possible. Significant dif- 
ferences between fit and data are only observed outside 
the signal region. The mass fit has 2905 ±58 A" — > A+7T - 
signal events, 252±46 other fully reconstructed A° candi- 
dates (which are also used to determine the A£ lifetime) , 
and 11% background in the signal region. 




5.2 5.4 5.6 5.8 6.0 6.2 

m(Aj it") GeV/c 2 



FIG. 1: The distribution of the invariant mass of A b — > At it 
candidates (points) with the fit overlaid (solid black line). 

Due to the trigger requirements on the track do and 
track-pair L xy , the observed A° ctdistribution is not a 
simple exponential. Consequently, an efficiency (e(ct)) 
must be included to model the acceptance of the trigger 
and offline selection. The largest corrections are due to 
the do requirements of the TTT. The two-dimensional 
ct — <j ct probability density function (pdf) for the signal 
and other fully reconstructed A° components is given by 

P(ct,cr,a ct ;S lt2 ) = P(ct\er,a ct , S lt2 )-P((Tct)-e(ct)- (1) 

where Si^ are the two descale factors obtained from a 
two-Gaussian modeling of the resolution function in MC 
(one for each Gaussian). The scale factor is necessary be- 
cause the kinematic fitter underestimates the uncertainty 
on the ct (the same scale factor is used for signal and back- 
ground). P(ct\a ct , £1,2) is a one-dimensional conditional 
pdf for observing this value of ctgiven the true A^ lifetime 
(t), a ct , and S\$. For the fully reconstructed A° com- 
ponents this pdf is a decreasing exponential convoluted 
with the sum of the two resolution Gaussians. P(cr ct ) is 
the pdf for observing a ct and is obtained from the side- 
band subtracted data distribution, where the sideband is 
defined as 5.8 < m(Ag) < 7.0 GeV/c 2 . For each back- 
ground component Eq. Q] is modified in a suitable way, 
apart from the partially reconstructed B mesons, which 
do not populate the signal region and are therefore not 
included in the lifetime fit. 

A sample of simulated signal events is used to extract 
e(ct). This sample consists of single b hadrons gener- 
ated with a pt spectrum extracted from the data sam- 
ple and decayed with EvtGen This MC sample is 
further reweighted in order to match the data in a num- 
ber of relevant variables: the choice of 'trigger tracks' 
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(the pair of final state particles which cause the TTT 
to fire), the proton production angle in rest frame 
which is sensitive to A° polarization, and the contribu- 
tions of the A+ Dalitz components 19]. The TTT effi- 
ciency function is represented by a histogram calculated 

as e(ct) = h(ct)/J2 l exp{ct,cT MC )®R(S 1:2 ,crlt)- Thc 
numerator is a smoothed histogram of the cifor all MC 
events that pass the trigger and analysis selection crite- 
ria. Each bin of the denominator is calculated by sum- 
ming the analytical cidistribution at the cibin center over 
all events (indexed by i) that pass the criteria required 
to fill the numerator. The analytical ctdistribution is an 
exponential convoluted with the resolution function R. 
Fig. [3 shows the resulting finely binned TTT efficiency 
histogram used for the A° signal components. Exactly 
the same procedure was used to derive an efficiency his- 
togram for the fully reconstructed B meson background. 




-0.05 0.00 0.05 0.10 0.15 0.20 0.25 0.30 

ct(A°) [cm] 



FIG. 2: The smoothed A° b TTT efficiency histogram (line) 
superimposed on the unsmoothed histogram (points) 

Our approach assumes that the simulation of trigger 
and detector can be used to derive e(ct). This assump- 
tion can be validated in data using J/tp — > fi + /i~ decays 
collected by the dimuon trigger which does not bias their 
lifetime. The observed four momenta of J/tp — > /i + [iT de- 
cays were also used as the input for a simulated sample 
of J/tp — ► /i + /i~ decays, subsequently fed to the TTT 
and detector simulation (data-seeded MC). Comparing 
the number of real J/tp — > /i + /i~ decays that pass the 
TTT with the number of data-seeded simulation decays 
that pass the TTT simulation gives a direct check of the 
reliability of the simulation. For both real and simulated 
J/tp decays we compute the TTT efficiency as a function 
of L xy and form their ratio R t {L xy ). The deviation of 
the slope of R e (L xy ) from is a measure of the quality of 
modeling of TTT in the simulation. Observed R e (L xy ) is 
incompatible with a null slope at the 3— 4<r level; we treat 
this discrepancy as a source of systematic uncertainty. 

We perform an unbinned maximum likelihood fit for 
cr(Ag) to the data that yields cr(Ag) = 420.1 ± 13.7 /on 
(stat). The total likelihood is: C = J], V, Nf g ■ 
Pj(cti, a c t i] 5*1,2), where the subscript i runs over events, 



and the subscript j runs over classes of event: the 
fully reconstructed signal and background fit compo- 
nents, the combinatorial background, and the partially 
reconstructed A° decays. Pj(cti, a c t %; £1,2) is a two- 
dimensional pdf of the form given in Eq. [TJ 7V| lg is thc 
number of events of this class occurring in the signal re- 
gion. The resulting likelihood projected onto the ci-axis 
is shown in Fig. 03 We fit the data for the A" lifetime 
after all procedures are established. The fit probability 
is estimated to be 37% using an unbinned Kolmogorov- 
Smirnov test. 




FIG. 3: The distribution of the ctof A b — > A^T -k~ candidates 
(points) with the fit projection overlaid (solid black line). The 
partially reconstructed A° fit components only have a 1% con- 
tribution to the signal region and are not shown. 

For each source of systematic uncertainty, we generate 
sets of events for about 500 pseudo experiments from a 
modified pdf and fit with both the standard and modi- 
fied fits. The mean of the distribution of thc difference 
between fit results obtained with the standard and mod- 
ified pdf is used as the systematic uncertainty. We con- 
sider the systematic uncertainties in two groups based on 
whether they affect the TTT efficiency or not. 

In the first group, the systematic uncertainty due to 
the alignment of the silicon detector is quoted from a 
previous study [5j (2.0 /im) where internal silicon sensor 
deformations and global misalignments of the silicon de- 
tector relative to the outer tracking volume are taken into 
account. The uncertainty due to the background compo- 
nent normalizations was taken into account by varying 
them according to their uncertainties derived from the 
mass fit (1.0 /im). 

In the second group of uncertainties, where the TTT 
efficiency is directly affected, the leading source is due 
to the slope of R e (L xy ) (8.6 /im). The uncertainty 
due to the A+ Dalitz structure is evaluated by vary- 
ing the relative contributions of each Dalitz compo- 
nent according to the world average uncertainty [161 ] 
(3.7 /im). The effect of an uncertainty in the combina- 
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torial background citemplate is computed by modifying 
it to a smoothed version of the actual upper sideband 
cidistribution (2.9 fim). The uncertainty due to the par- 
ticle identity of the tracks which fired the trigger is eval- 
uated by varying the relative contributions of different 
trigger-track combinations in the MC (2.0 fim). The un- 
certainty due to the A° polarization is obtained by vary- 
ing the slope of the MC reweighting factor by one-sigma 
from a straight line fit for the proton production angle in 
the A+ rest frame (1.4 /jm). The uncertainty due to the 
transverse position of the pp primary interaction point is 
computed by dividing the MC into independent subsam- 
ples representing the extreme variations of the primary 
interaction point (1.2 fim). The uncertainty due to the 
TTT efficiency used for the B° background is evaluated 
by tightening the mass cut on the D + tt~ candidate in 
the underlying B° MC reconstruction (1.0 fim). The un- 
certainty due to the lifetime assumed for the B° back- 
ground is obtained by varying this lifetime according to 
the world average uncertainty [1 61 ] (1.0 /mi). The uncer- 
tainty due to a correlation between the do requirements in 
the SVT and reconstruction levels is estimated by smear- 
ing the latter in the signal MC by an amount extracted 
by comparing their difference distributions between data 
and MC (1.0 /urn). The total systematic uncertainty is 
computed by adding all the contributions in quadrature, 
which is 10.6 fim. 

Numerous cross-checks were performed. We used our 
procedure to measure the B° lifetime in the B° — > 



D 



and B° -> D~ 



decay modes and the B + life- 
time in the B + — > D 7r + mode. These B° and B + life- 
time measurements are statistically consistent with the 
world averages [l6[- We checked the effect of uncertain- 
ties in the mass template shapes, the pt spectrum of the 
A° , the effect of assuming different A° lifetimes in the 
MC, the scale factor applied to the cr ct , the A+ lifetime, 
and the model of the uncertainty of the transverse posi- 
tion of the pp primary interaction point. We also checked 
the effect of the uncertainty in the shape of cr ct and used 
a large signal MC sample to verify that the fitter itself 
does not introduce a bias in the measured lifetime. 

In summary, using a sample of 2905 ± 58 fully recon- 
structed A" — > A+7T - decays we measure the lifetime 
of the A° baryon to be t(A°) = 1.401 ± 0.046 (stat) ± 
0.035 (syst) ps (corresponding to cr(A") = 420.1 ± 
13.7 (stat) ±10.6 (syst) fim where c is the speed of light). 
This is the single most precise measurement of the A° 
lifetime. 

Using the current world average for the B° life- 
time 0, we obtain r(Ag)/r(B ) = 0.918 ± 0.038 (stat 
+ syst). There is good agreement between our re- 
sult and the current world average of r(A°)/r(i? ) = 
0.99 ±0.10 16], and between our result and the previous 
CDF result [5j. This measurement is also compatible with 
the current HQE value [3] of r(A°)/r( J B ) = 0.86 ±0.05, 
thus supporting the HQE picture of weak decays of heavy 



baryons. 
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